Performance of solenoids vs. 



quadrupoles in focusing and energy selection of laser 
accelerated protons. 
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Using laser accelerated protons or ions for various applications - for example in particle therapy or 
short-pulse radiographic diagnostics - requires an effective method of focusing and energy selection. 
We derive an analytical scaling for the performance of a solenoid compared with a doublet /triplet 
as function of the energy, which is confirmed by TRACEWIN simulations. The scaling shows that 
above a few MeV a solenoid needs to be pulsed or super-conducting, whereas the quadrupoles can 
remain conventional. The transmission of the triplet is found only 25% lower than that of the 
equivalent solenoid. Both systems are equally suitable for energy selection based on their chromatic 
effect as is shown using an initial distribution following the RPA simulation model by Yan et al.Q. 



I. INTRODUCTION 

Laser acceleration of protons or ions requires ultra 
high laser intensities focused on thin target foils. This 
has been demonstrated in numerous experiments (see, 
for instance, Refs. [2|-[5j). Applications of this novel ac- 
celeration method have been suggested, for example in 
terms of a new proton source for radiation therapy 
as alternative to more conventional accelerator technolo- 
gies, like protons or ions from cyclotrons or synchrotrons. 
Other potential applications might be proton radiogra- 
phy, neutron imaging or isotope production. Energies up 
to 170 MeV for deuterons have recently been observed at 
the TRIDENT laser and explained as "break-out after- 
burner" (BOA) mechanism [9], with possible applications 
as very short-pulse neutron source. 

All of these applications have to cope with the charac- 
teristics of laser accelerated protons or ions: a large en- 
ergy spread as well as angular spread, sub-ns time scales, 
significant shot-to-shot fluctuations and - for practical 
applications - relatively low repetition rates compared 
with conventional accelerators. This requires specific 
methods to suitably manipulate laser accelerated protons 
in space and time and match them to the need of an ex- 
periment or application. 

In a preceding study we have shown that a single 
solenoid magnet can be used very effectively to combine 
angular focusing (collection) with energy selection due to 
the lens chromatic effect [lOj. This is extended here to a 
comparative evaluation of quadrupole focusing (doublet 
or triplet) and solenoid focusing. As in the solenoid case, 
we use the dependence of focal length on energy and em- 
ploy a radially confining aperture to select a suitable en- 
ergy window. The main difference of quadrupoles versus 
a solenoid is their first order focusing property (solenoids 
focus in second order) and the asymmetry in focusing 
(astigmatism) as well as different chromatic effects be- 
tween the horizontal and vertical planes. On the other 



hand alternating gradient focusing with quadrupoles is 
known to be more efficient at increasing energy and helps 
to avoid pulsed magnet or superconducting technology. 

It should be mentioned here that the alternative to 
energy selection by the chromatic effect is the more con- 
ventional energy selection using the dispersive properties 
of dipole magnets. All therapy oriented studies on laser 
acceleration carried out in the course of the last years 
have used such a system with a collimating aperture fol- 
lowed by a dispersive dipole magnet [lll-KUj. Dispersive 
energy selection is certainly an option, but schemes with- 
out transverse focusing may result in significant efficiency 
loss. Chromatic energy selection, instead, combines the 
advantages of focusing with energy selection in a single 
device. However, for effective selection it is required that 
the beam size by chromaticity dominates over the size 
by the intrinsic emittance, which is not the case in many 
situations - but certainly for laser accelerated particles. 

An interesting alternative is the micro- lens, which fo- 
cusses by means of an electron plasma initialized by a sec- 
ond short-pulse laser beam and collapsing inside a cylin- 
der [15j| . Energy selection is enabled by using time-of- 
flight and triggering the second laser beam accordingly. 
For therapy applications, however, where energy selec- 
tion failure cannot be tolerated, the fully controllable 
powering of conventional magnets for energy selection 
seems essential. 

In Section [II] we compare solenoidal with dou- 
blet/triplet focusing. Their chromatic focusing proper- 
ties are studied in Section HIB In Section HVl results are 
applied to energy selection by an aperture for both sys- 
tems using an input distribution from the simulation of 
a radiation pressure acceleration model. In Section IVl we 
draw conclusions. 



II. COMPARISON OF SOLENOIDAL AND 
QUADRUPOLAR FOCUSING 
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Solenoids are frequently applied in injectors for fo- 
cusing of particles with relatively low energy and large 
divergence. Likewise, focusing of laser accelerated pro- 
tons with energies of the order of 10 MeV was demon- 
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strated using a pulsed solenoid [H, [l7| . Other laser pro- 
ton experiments in a comparable energy range have suc- 
cessfully employed small aperture, high-gradient perma- 
nent magnet quadrupoles [18[ . Preference of quadrupoles 
over solenoids depends on the individual application, but 
certainly energy and the question of room-temperature, 
non-pulsed quadrupoles versus pulsed or superconduct- 
ing solenoids matter. 

A useful guidance to decide on the basis of required 
field strengths can be obtained from a scaling expression 
using the thin lens approximations for the focal length 
f s of a solenoid and Fd of a quadrupole doublet as sug- 
gested in Ref. [19j. With B the field strength (for the 
quadrupoles defined at the poles), and assuming that 
both focal lengths are defined from the respective cen- 
ters, we have for the solenoid of length L 
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where I is the individual quadrupole length, s the sepa- 
ration of quadrupoles (from center to center) and a the 
maximum beam radius (pole radius). Note that the fo- 
cal strength of a doublet increases with the separation of 
its components - of course on the expense of decreasing 
acceptance. Comparing a solenoid with a doublet of the 
same overall length L and equal field B, we readily ob- 
tain from Eqs. [T] and [2] the ratio Td of focusing strengths 
(here defined as inverse focal lengths) in terms of only 
geometrical quantities: 
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Eq. [3] indicates that the focusing strength of a doublet 
is superior to that of a solenoid, if a is sufficiently small 
relative to the length. As an example, consider a doublet 
with a gap between magnets equal to their length, in 
which case we have Td = (2/3) 3 (L/a) 2 and the transition 
condition T d > 1 occurs for L/a > (3/2) 3/2 . 

In order to obtain systems with equivalent focal lengths 
it may be required to adjust B s according to Eq. [TJ This 
results in an effective field for the solenoid 
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which may lead to the requirement of super-conducting 
or pulsed power technologies for the solenoid. 

We can apply this to the collection of laser particles, if 
we assume the focal spot is at the source and the beam 
is to be made parallel by the lens. In the interest of 
smoother focusing we find it preferable to use a triplet 
rather than a doublet as reference case. In Fig.[T]we show 
- as example - a triplet with L = 0.334 m, I = 0.06 m 
and gradients of 30 T/m, -30 T/m and 15 T/m. The 



calculation of matched beam optics is obtained using the 
envelope option of the TRACEWIN code [20], which is 
also employed further below for particle tracking. 

2 MeV protons with source divergence of ±125 mrad 
are made parallel with a triplet focal length (source to 
triplet center) F t =0.254 m, maximum envelope a = 48 
mm, hence maximum pole-tip field B t = 1.44 T. Note 
that the distance source to lens is given by F t — L/2, which 
is 87 mm in this case. The energy spread is assumed to be 
zero, which is to avoid the chromatic energy effect at this 
point. For comparison we also show a solenoid focusing 
with the same L; in order to achieve the same focal length 
we require a slightly increased magnetic field B s = 1.53 
T. Next we check if our TRACEWIN results still obey the 




FIG. 1. TRACEWIN envelopes for equivalent triplet (top) 
and solenoid (bottom) solutions at 2 MeV. 

scaling of Eq. [3j although we have replaced the doublet 
by a triplet. We find good agreement, if we replace 4s by 
2s. The need for this adjustment can be interpreted as 
necessary compensation of the (approximate) doubling 
of the triplet length compared with that of the doublet. 
This suggests a triplet focusing enhancement factor 
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Assuming that the gap equals the quadrupole length, 
we have 5 = 2/, which results in T t « 1.12 in the above 
example. Applying Eq. [4] we find that the predicted 
solenoid field for equivalent focal length is 1.52 T, which 
agrees quite well with the above TRACEWIN result - 
in spite of the thin lens approximations employed in the 
derivation of T t . 

In order to further examine the validity of Eq. [5] for 
different energies we extend the systems in Fig. [1] to pro- 
ton energies of 0.2, 20 and 200 MeV, again assuming 
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at each energy equal overall lengths L for solenoid and 
triplet, equal focal lengths for the two systems as well 
as equal gap and quadrupole lengths in the triplet case. 
For the initial divergence we assume - somewhat arbi- 
trarily - a divergence scaling x' oc (/#7) -1 / 2 to account 
for the expected trend of decreasing divergence with en- 
ergy. The value of x' at 2 MeV is kept as before. We first 
use TRACEWIN and search again for matched solutions 
requiring a parallel output beam for vanishing energy 
spread. In Table [IT] we summarize all relevant parame- 
ters including the resulting quadrupole pole tip fields B t 
and solenoid strengths B s . The focal length F is again 
defined from source to center of the respective lens sys- 
tem and found to increase with energy. The theoretically 
expected triplet enhancement factor T t is calculated from 
Eq. [5] (s = 21) by inserting the respective geometrical di- 
mensions. Using Eq. [4j which applies equally to the pole 
tip field of the triplet as it does for the doublet, we can 
thus derive the theoretically expected B s * and compare 
it with the actual B s obtained from TRACEWIN match- 
ing. As result we find an overall good agreement between 
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TABLE I. Comparison simulation - theory for equivalent 
solenoid and triplet focusing properties; lengths in cm and 
magnetic fields in T (x' in mrad). 



the TRACEWIN calculated B s and the theoretical B*, 
which confirms the theoretically derived triplet focusing 
enhancement over a solenoid. In Fig. [2] we summarize 
the main findings from this comparison. 
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FIG. 2. Comparison of triplet and solenoid systems with equal 
focal lengths as function of energy. 

In summary, this demonstrates that for sub-MeV or 
few MeV energies solenoids are a convenient approach, 
whereas the quadrupole doublet /triplet (or multiplet) 
has advantages for higher energies as its pole tip field 
strengths remain within iron saturation. Eq. [5] also sug- 



gests that T t oc 1/a 2 , hence larger a due to an in- 
creased initial divergence shifts the advantage of the dou- 
blet / triplet over the solenoid to higher energies - and vice 
versa. 



III. CHROMATIC PROPERTIES OF SOLENOID 
AND TRIPLET 

TRACEWIN is used here for particle tracking. Al- 
though primarily a linear accelerator design and verifi- 
cation tool, it has a number of features, which make the 
code suitable for our problem as well. In particular, it 

• is a self-consistent 3D particle-in-cell code suitable 
for tracking of up to 10 T simulation particles, 

• is capable of energy dependent focussing (chro- 
matic aberrations), 

• includes field map options for magnetic elements 
(like solenoids) to model higher order (for example 
geometric) lens aberrations, 

• provides standard 6D phase space initial distribu- 
tions as well as user provided input distributions; 
as "standard uniform" we choose for this study the 
option of a uniform distribution in the 4D trans- 
verse space space as well as uniform within the lon- 
gitudinal phase plane ellipse; 

• includes an envelope option for beam optics design. 

Space charge options with 2D/3D Poisson solvers exist, 
but space charge is ignored in this study. In Ref. [lOj it is 
shown for solenoids that space charge is generally weak; 
in the near-source region, where extremely high proton 
densities are prevailing, neutralization by the co-moving 
electrons helps. 

In the following we use as reference a solenoid and an 
equivalent triplet, both 104 cm long and designed to bring 
250 MeV protons to a focus at 2.73 m with the following 
assumptions: 

• initial maximum divergence angle: ±28 mrad 

• energy spread: practically mono-energetic 

• distance laser target - first magnet: 35 cm 

• beam pipe radius: 3.5 cm 

• aperture radius of solenoids and quadrupoles: 3 cm 

• length of solenoid field map: 104 cm 

• length of solenoid field region: 80 cm 

• averaged solenoid field 6.27 T 

• length of solenoid field map: 96 cm 

• length of quadrupoles: 15 cm 
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• quadrupole pole tip fields: 1.5/1.5/1.0 T 

Fig. [3] shows density plots from a multi-particle simu- 
lation using a low number of simulation particles (only 
3000), which helps to visualize single particle rays. The 
simulation was carried out with the "standard uniform" 
initial distribution of TRACE WIN. Maximum energy de- 
viations in this example have been chosen as ±5 x 10 -4 
MeV centered at 250 MeV, hence practically mono- 
energetic. The common waist for x and y for the triplet 
(stigmatic image) is relevant for optimum energy selec- 
tion as will be shown in the next section. It is noticed 
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FIG. 3. Density plots for TRACEWIN obtained reference 
cases for "equivalent" solenoid (top) and triplet (center:x and 
bottom:y) focusing. 



to the geometrical mean of the two coefficients for the 
equivalent triplet. 

The chromatic effect is strongly correlated with time 
as higher energy particles travel ahead. The stronger fo- 
cusing for lower energies leads to an enhanced transverse 
phase space rotation of protons at the bunch end as com- 
pared with the high energy particles at the bunch head. 
The resulting slip in the transverse phase planes causes 
an effective transverse emittance increase, which can sig- 
nificantly exceed the initially small production emittance. 
The effective emittance is obtained by averaging the in- 
stantaneous emittances over the full bunch length, hence 
the full energy spectrum fioj ] . 

IV. TRANSMISSION AND ENERGY 
SELECTION 

As suggested in Refs. [10| the pronounced chromatic fo- 
cusing effect can be used for an effective energy selection, 
if the beam is focused into a suitably defined transverse 
aperture. Only particles with focal spot sufficiently close 
to the aperture plane are transmitted effectively. As for 
solenoids this works effectively only if the beam is "chro- 
maticity dominated": at a selection aperture the beam 
size by chromaticity dominates over the size generated by 
the intrinsic emittance at any relevant value of the en- 
ergy (i.e. position along the bunch). This is always the 
case for laser generated ions with their extremely small 
emittance at any given energy, which is owed to the very 
small source spot size. It should be mentioned here that 
the intrinsic emittance should also include emittance in- 
crease due to higher order aberrations of the lens, which 
can be a problem for short solenoid lenses at low energy 
but is not further considered in this study. 

A. RPA generated initial distribution 



that the de-focusing effect of the first quadrupole in y 
leads to a - in this example - small beam loss at the 
aperture of the second quadrupole. 

In Refs. [l0[ the energy dependence of the focal length 
of a solenoid lens was expressed in terms of a chromatic 
coefficient, which we generalize here to cope with the 
different focusing in x and y for a triplet: 



_ fifx,y/ fx,y 

5E/E ' 



(6) 



Here f XiV is the focal length at the reference energy E 
and a x ^ y is specific to the geometry of the focusing set- 
up. For the examples of Fig. Owe find from TRACEWIN 
simulation for the solenoid a r « 1.9 and for the triplet 
a x « 0.9 as well as a y « 3.8. The much larger a y is 
a result of the de-focusing in y at the first lens and the 
thus much larger overall envelope excursions in y. The 
solenoid chromatic coefficient is - not surprisingly - close 



For practical considerations it is advantageous to con- 
sider an initial distribution in 6D phase space with a 
broadened energy distribution according to some laser 
acceleration model. The Radiation Pressure Acceleration 
(RPA) mechanism [21-26] has a high potential to reach 
proton energies of hundreds of MeV. A specific theoreti- 
cal version of it has been discussed in Ref. UJ and applied 
to proton therapy conditions in Refs. [To|, H?) to create a 
proton energy spectrum extending up to 250 MeV. The 
output of this RPA-simulation can be described as spec- 
tral yield 



dN(E, g) 
dE 



[MeV~ 



(7) 



which describes the number of particles in an energy in- 
terval dE and within a cone angle ±fi. The thus defined 
proton spectrum (details see Ref. [27[) is plotted in Fig.HJ 
Its energy distribution is peaked above 200 MeV - deter- 
mined by the laser intensity - with a relatively broad foot 
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towards lower energies. As input into our TRACE WIN 
simulations we take a bi- Gaussian approximation to this 
energy spectrum shown by the continuous curve in Fig.SJ 
The 6D initial distribution is taken as a Gaussian ran- 
dom distribution in the variables t,x,x' ,y,y' . For the 
rms widths in x' , y' we have chosen 35 mrad - in contrast 
with the broader tails in divergence indicated in Fig. [4j 
which are probably due to the 2D nature of the RPA sim- 
ulation. The initial spot radius and pulse duration are in 
the /am rsp. ps scales; their actual values play no role as 
long as space charge is considered as neutralized initially. 
We also note that the detailed profile of the energy spec- 
trum is only exemplary - what matters primarily is its 
gradient near a selected energy. 



2.5 



.x 10 



Z 

T3 



- 10mrad 
-20mrad 
-30mrad 

40mrad 
-50mrad 
-75mrad 

100mrad 





100 150 200 

proton energy[MeV] 



250 



FIG. 5. Density plots for RPA-distribution in equivalent 
solenoid (top) and triplet (centenx and bottom:y) systems 
adjusted to 220 MeV. 



FIG. 4. Spectral yield of protons as a function of energy 
and for different capture cone angles Q, with bi- Gaussian fit 
(continuous line). 



B. Comparative transmission 

For this purpose we reduce the magnet fields in the 
equivalent solenoid and triplet systems of Fig. [3] for nom- 
inal transmission at 220 MeV, which is closer to the peak 
of the energy spectrum. We also require - arbitrarily - 
a focus at the distance of 2.73 m from the laser target, 
where the energy selection aperture is placed. Employ- 
ing the above defined RPA-distribution and a 3 mm ra- 
dius aperture, the resulting orbits of a TRACEWIN sim- 
ulation with 3000 rays are shown in Fig. In Fig. [6] 
we examine the transverse emittances for the triplet of 
Fig. The relatively large spread in x' \y' and energy 
width together with the energy dependent focusing re- 
sult in significant emittance growth in x and y within 
the quadrupoles, accompanied by emittance reductions 
due to beam loss on the radial aperture. Note that 
emittances are understood here as averaged over the full 
bunch length, while the "instantaneous" emittances at a 
given position along the expanding bunch remain at their 
small initial values. The emittance reduction in y within 



the first quadrupole reflects the beam loss in the defocus- 
ing y— direction. A comparison of the equivalent solenoid 



[1/8/2013] TraceWin - CEA/D5M/Irfu/SACM 




Position (m) 

FIG. 6. Transverse (normalized and bunch averaged) rms 
emittances for triplet case. 

and triplet focusing systems shows that the overall trans- 
mission for the solenoid is 47%, and 35% for the triplet. 
The corresponding loss profiles are shown in Fig. [71 The 
beam loss in y in the first lens of the triplet is to some 
extent compensated by an enhanced transmission in x, 
which explains why the triplet transmission is only 25% 
lower than the solenoid one. A more differentiated in- 
sight is gained, if we truncate the divergence of the initial 
Gaussian distribution in x' ', y' by eliminating all particles 



6 



80 




Position ( m ) 




Position ( m ) 

FIG. 7. Loss profiles for RPA-distributions in the equivalent 
solenoid (top) and triplet (bottom). 



beyond a "divergence limit" as shown in Fig. [51 Below 
about 40 mrad the solenoid accepts all injected particles 
in the truncated distribution (intensity in it relative to 
un-truncated distribution indicated by the dotted line); 
at this value the acceptance limit is reached and larger 
divergence particles are lost. The triplet, instead, starts 
losing particles above 20 mrad, but the transition to its 
acceptance limit is smoother - apparently due to the ben- 
efit from the horizontal plane. 
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FIG. 8. Transmission of solenoid and triplet as function of 
an upper cut-off (divergence limit) of the injected particles. 
The transmission is in % of the particles in the original un- 
truncated distribution. 



and chromatic coefficient a, 

Ra = a—^-Am ax , (8) 

where A max is the maximum envelope at the lens. Using 
Amax ~ 3 cm and a^2we expect that an energy width 
of ±4% (±8.8 MeV) should be obtainable with an aper- 
ture of 2.4 mm radius. This is approximately confirmed 
by the energy spectra in Fig. [9j where Ra was chosen 
as 3 mm for the solenoid and 2.7 mm for the triplet to 
reach the same fwhm width. The overall yield in the 
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FIG. 9. Selected energy spectra for equivalent solenoid (Ra = 
3 mm) and triplet systems (Ra — 2.7 mm). Top: initial; 
center: solenoid; bottom: triplet. 

selected energy windows is 17% for the solenoid and 13 
% for the triplet, which follows approximately the 25% 
triplet transmission reduction found in Fig. [51 

For the triplet we have also simulated an elliptical se- 
lection aperture in x, y with the same area but semi-axes 
in the ratio 1:4 to match the ratio of chromatic coeffi- 
cients according to Eq. [8j However, we obtain practi- 
cally the same energy selection width and profile as well 
as transmission. This seems somewhat unexpected, but 
apparently the loss of selection in one plane is compen- 
sated by better selection in the other plane. 



C. Selection of energies 



V. CONCLUSION 



Following Ref. [10| the radius of a selection aperture 
is proportional to the product of required energy width 



The purpose of this study has been a comparative as- 
sessment of the focusing properties of a solenoid and a 
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quadrupole triplet in the context of laser accelerated pro- 
tons (or ions). Possible application of this acceleration 
method can be envisioned in the field of particle therapy; 
but also in proton radiography in areas, where energies of 
a few hundred MeV, moderate integrated intensities but 
high peak intensities in time-scales of few ns or sub-ns are 
needed. The relatively large initial angular and energy 
spreads are a challenge for all of these applications. 

In terms of transmission it is found that "equivalent" 
systems - same geometrical length and apertures - give 
only the relatively small reduction in transmission of 25% 
for the triplet vs. the solenoid, which is owed to the 
un-symmetric focusing of quadrupoles. For increasing 
energies - already above a few MeV - the weaker fo- 
cusing properties of solenoids require pulsed or super- 
conducting technology, whereas pole-tip fields of a dou- 
blet or triplet can remain within room-temperature iron 
saturation. This appears to be a clear advantage for 
quadrupoles in future therapy applications, where short- 
term and fully controlled changes of energy and magnetic 



rigidity are required. 

The large energy spreads lead to a dominance of chro- 
matic effects, which can be used for energy selection. 
Equivalent solenoid and triplet systems are equally suit- 
able for this purpose in spite of the strongly differing 
chromatic coefficients in x and y. Chromatic effects lead 
to inevitable correlations between energy and transverse 
position, which cannot be ignored for therapy applica- 
tions. In Ref. [10] it is shown that properly placed scatter 
targets can be used to remove these correlations. 

The role of space charge and geometric aberrations - 
dominant for short solenoids - needs further consideration 
even though they are not expected to alter the major 
conclusions. Noting that solenoid focusing is independent 
of the charge, the neutralizing co- moving electrons, which 
are always present in laser acceleration, will be strongly 
focussed towards the axis [28|. In the triplet, instead, 
these electrons will be de-focussed to the aperture when 
entering the first quadrupole, which may have an effect on 
the quality of focusing and needs to be further explored. 
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